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Abstract: A wide range of atmospheric monitoring instruments is employed at the Pierre Auger Observa- 
tory : two laser facilities, elastic lidar stations, aerosol phase function monitors, a horizontal attenuation 
monitor, star monitors, weather stations, and balloon soundings. We describe the impact of analyzed 
atmospheric data on the accuracy of shower reconstructions, and in particular study the effect of the data 
on the shower energy and the depth of shower maximum (Xmax)- These effects have been studied using 
the subset of "golden hybrid" events — events observed with high quality in the fluorescence and surface 
detector — used in the calibration of the surface detector energy spectrum. 



Introduction 

The Pierre Auger Observatory in Malargtie, Ar- 
gentina is a hybrid facility that uses its atmospheric 
fluorescence detector (FD) to obtain calorimetric 
estimates of shower energies. The atmosphere acts 
both as a calorimeter and a scintillator. It affects 
the fluorescence yield of showers and attenuates 
Ught between showers and the FD. The atmosphere 
over the observatory is also dynamic. Therefore, 
the state of the atmosphere must be continuously 
monitored to ensure reliable energy estimates. 

To an excellent approximation, the molecular and 
aerosol scattering processes that contribute to the 
overall attenuation and scattering of light in the at- 
mosphere can be treated separately. In Malargtie, 
the molecular component is determined by regular 
measurements of several macroscopic parameters, 
including altitude profiles of air temperature, pres- 
sure, and density. It has been shown that daily vari- 
ations in these parameters have a small impact on 
shower energy estimates (AE/E < 1%) and the 
depth of shower maximum (AX^ax — 6 g cm~^); 
hence, the observations have been incorporated 
into monthly models for use in the FD reconstruc- 
tion [1]. A more important factor is the effect of 
humidity on the fluorescence yield, because the 
yield provides a scaling factor for the energy. Pre- 



liminary estimates suggest a 5 — 10% effect on the 
fluorescence yield near the ground, and < 3% at 
4 km above sea level [2]. This effect has been in- 
corporated into the reported uncertainty of the flu- 
orescence yield [8]. Aerosols, unlike the molecu- 
lar component of the atmosphere, are much more 
variable, and can change significantly in the course 
of a few hours. Therefore, aerosols are systemat- 
ically measured at all FD sites, and the parame- 
ters most important for the FD reconstruction are 
recorded hourly. The aerosol data and the facilities 
used to collect them are described in detail else- 
where [3]. In this paper, we discuss the effect of 
aerosol measurements on energy and Xj^ax for an 
important subset of observed showers. 

Aerosol Measurements 

The presence of aerosols does not influence the 
air fluorescence yield, so their primary effect on 
the shower reconstruction comes from their role in 
Ught attenuation and scattering. Fluorescence light 
is emitted isotropically, so a detector with a field 
of view AQ observing shower light of intensity Iq 
will observe a light level 

AQ 

I=Io-Tm-Ta-{l + H.O.) ■ — (1) 
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In this expression, Tm is the transmission factor 
due to molecular scattering; Ta is the transmission 
due to aerosol scattering; and H.O. is a higher- 
order correction that accounts for the single and 
multiple scattering of photons into (or out of) the 
detector field of view. Ignoring multiple scattering 
effects, Ta and H.O. may be fully characterized 
using three independent measurements: the height 
profile of the vertical aerosol optical depth (VAOD) 
r(/i); the wavelength dependence of the VAOD; 
and the normalized aerosol differential scattering 
cross-section, or phase function, P{9). 

The VAOD primarily affects light attenuation. It 
is defined for each height above the ground level, 
so that aerosol transmission Ta{h) between the 
ground and height h is 

Ta{h) = e-^(''). (2) 

The height profile r(/i) is measured independently 
by three elastic backscatter lidar stations and FD- 
reconstructed laser tracks provided by the Central 
Laser Facility (CLF) [3, 4, 5]. With a typical mea- 
surement uncertainty of ±0.01, the CLF measure- 
ments are currently used for shower reconstruc- 
tions. However, both the lidars and the CLF use 
monochromatic light sources, so the wavelength 
dependence of the VAOD is measured by two other 
independent instruments: the Horizontal Attenua- 
tion Monitor (HAM) and the robotic astronomical 
telescope FRAM [6, 7]. This wavelength depen- 
dence can be parametrized in terms of the so-called 
Angstrom exponent 7: 

r(A) = ro.(^)\ (3) 

where tq is measured at a reference wavelength 
Ao = 355 nm. Observations made between June 
and December 2006 indicate a typical value 7 = 
0.7 ± 0.5 for the observatory location. 

Finally, aerosols not only attenuate light from 
air showers, but also scatter Cherenkov light into 
the FD field of view, contaminating the fluores- 
cence signal. The angular distribution of aerosol- 
scattered light is given by the aerosol phase func- 
tion, which we model using two free parameters 
(see [3]). The parameter / is sensitive to the rel- 
ative strength of forward and backward scatter- 
ing and g, the mean cosine of the scattering an- 



gle, is the measure of the asymmetry of scattering. 
Aerosol Phase Function monitors (APFs) perform 
hourly measurements of these aerosol scattering 
properties at two FD locations [3]. 

Analysis 

The showers of greatest importance to Auger are 
those measured with high quality in both the FD 
and the surface detector (SD), because these are 
used to set the energy scale of the overall detector. 
For the analysis presented here, the set of events 
passing strong quality cuts presented in [8, 9] are 
used to determine the effect of aerosol measure- 
ments on energy and X,^^^ estimated by the FD re- 
construction. Only cloud-free measurements, iden- 
tified by a strict quality cut on the FD longitudinal 
profile, were used in this study. 

Several different studies are of interest: 

L The use of aerosols in the reconstruction, com- 
pared to the use of a pure molecular atmosphere. 

2. The propagation of measurement uncertainties 
in VAOD, 7, /, and g in the FD reconstruction, 
and in particular their effect on energy and X^ax- 

3. A test of the assumption of atmospheric hori- 
zontal uniformity used in the FD reconstruction. 

(1) Effects of the Presence of Aerosols on Air 
Shower Detection 

We have compared the reconstruction of show- 
ers using real-time aerosol measurements with the 
same events reconstructed using a purely molecu- 
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Figure 1: Energy difference for events recon- 
structed with measured aerosol parameters (Eaero) 
and with assumption of no aerosols (molecular at- 
mosophere; Emoi). 
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Figure 2: Spread in E and X^ax due to statistical fluctuations in aerosol measurements. AE_|_ and AX+ are 
calculated using the +1 a uncertainty (on aerosol measurement), AE_ and AX^ using the — la uncertainty. 



lar atmosphere. As shown in Figure 1, neglecting 
the presence of aerosols causes, on average, a 12% 
underestimate in shower energies. Moreover, the 
long tail in the distribution indicates the enormous 
effect of aerosol attenuation on the reconstruction 
of a significant fraction of all showers: 15% of 
the showers have an energy correction greater than 
25%; 6.5% of events more than 30%; and 3% more 
than 40%. Especially for the highest energy events, 
where the statistics are poor, real-time atmospheric 
calibration is essential. 

Having established the significant influence of 
aerosols on the reconstruction, we can also ask 
if simple parametric models of the aerosol con- 
tent are sufficiently accurate. For example, a 
two-parameter exponential aerosol density profile, 
characterizing conditions at the site, has been con- 
sidered. Preliminary studies indicate that this 
parameterization leads to a 4% overestimate in 
shower energies (with a large spread of 10%) 
compared to reconstructions performed with true 
aerosol measurements. 

(2) Uncertainties Introduced by Aerosol Mea- 
surements 

We have propagated the measurement uncertain- 
ties in the VAOD, phase function, and wavelength 
dependence in the hybrid reconstruction. Figure 
2 depicts the contribution of each measurement to 
the uncertainty in energy and shower maximum. 



The VAOD provides the dominant contribution to 
the shower uncertainties: 5.5% for the statistical 
uncertainty in energy and 4 g cm^^ for X^ax- The 
wavelength dependence and the phase function are 
significantly less important, contributing 1% and 
1.3%, respectively, to the uncertainty of the energy, 
and ^ 2 g cm"'^ to the uncertainty in X„ax- 
(3) Evaluation of the Horizontal Uniformity of 
the Atmosphere 

The atmospheric measurements at the Auger Ob- 
servatory, while extensive, are only able to ob- 
serve conditions at several locations across the site. 
Therefore, during the reconstruction of events one 
must assume that these limited measurements char- 
acterize broad regions of the atmosphere around 
each FD, or in other words, that the atmosphere 
exhibits a large degree of horizontal uniformity. 

The assumption of horizontal uniformity and its ef- 
fect on the reconstruction can be tested for both 
the molecular and aerosol components of the atmo- 
sphere. Molecular conditions are observed by bal- 
loon flights and two ground-based weather stations 
located at Los Leones and the Central Laser Facil- 
ity. Figure 3 indicates differences in the weather 
conditions observed at these sites, and the corre- 
sponding effect of these differences on shower re- 
constructions is < 1% for the shower energy and 
1 g cm"2 foj. x^^^. 
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Figure 3: Comparison of monthly atmospheric 
depth measurements at two sites separated by ~30 
km. 



For the aerosol component of the atmosphere, it is 
possible to estimate the effect of horizontal nonuni- 
formities by reconstructing showers using different 
sets of aerosol measurements. For this study, we 
have reconstructed a set of high quality air showers 
observed by the Coihueco FD using CLF VAOD 
profiles measured concurrently at Coihueco and 
Los Leones (the distance between these two FD 
sites is about 40 km). As shown in Figure 4, the 
systematic uncertainty on shower energies intro- 
duced by the assumption of uniformity is ~ 2.5%, 
with measurement fluctuations contributing ~ 7% 
to the statistical uncertainty. The shower maximum 
Xmax is also shifted by ^ 2.5 g cm^^, with typical 
statistical uncertainties of ~ 9 g cm~^. 

Discussion 

The following table summarizes atmosphere- 
induced uncertainties in the hybrid reconstruction: 



Effect 


AE/E 


A Xmax 


Molecular: 






Horizontal uniformity 


< 1 % 


1 gcm~^ 


Variations in air density profile 


1 % 


6 g cm~^ 


Aerosols: 






Horizontal uniformity | systematic] 


2.5 % 


2.5 g cm~^ 


Horizontal uniformity [statistical] 


7 % 


9 g cm~ 


Vertical Aerosol Optical Depth 


5.5 % 


4 g cm~^ 


Wavelength Dependence 


1 % 


1 gcm~^ 


Differential Scattering Cross-section 


1.3 % 


2 g cm~^ 
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Figure 4: Comparison of energy difference (top) 
and of shift in position of shower maximum (bot- 
tom) due to the use of non-local aerosol atmo- 
spheric parameters. 
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Table 1: Summary of atmosphere-induced uncer- 
tainties for the set of calibration events 



